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Parity-odd domains, corresponding to non-trivial topological solutions of the QCD vacuum, might 
be created during relativistic heavy-ion collisions. These domains arc predicted to lead to charge 
separation of quarks along the system's orbital momentum axis. We investigate a three particle 
azimuthal correlator which is a "P even observable, but directly sensitive to the charge separation 
effect. We report measurements of charged hadrons near center-of-mass rapidity with this observable 
in Au+Au and Cu+Cu collisions at ^/sNN =200 GeV using the STAR detector. A signal consistent 
with several expectations from the theory is detected. We discuss possible contributions from other 
effects that are not related to parity violation. 
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PACS numbers: ll.30.Er, 11.30.Qc, 25.75.Ld, 25.75.Nq 



Parity (V) violation in the weak interaction was ob- 
served in 1957 [1]. However, until recently, parity has 
been thought to be conserved in the strong interaction. 
Modern QCD theory does allow for parity violation, but 
experiments have not seen this violation and the result- 
ing constraints are tight [H [3] . Recently, it has been sug- 
gested that the hot and dense matter created in heavy-ion 
collisions may form metastable domains where the parity 
and time-reversal symmetries are locally violated j4|. In 
non-central collisions, these domains may manifest them- 
selves by giving positively and negatively charged parti- 
cles opposite-direction momentum 'kicks' along the an- 
gular momentum vector of the collision. The resulting 
charge separation is a consequence of two factors [IHl]: 
the difference in numbers of quarks with positive and neg- 
ative chiralities due to a non-zero topological charge of 
the metastable region, and the interaction of these par- 
ticles with the extremely strong magnetic field produced 
in such a collision (leading to the effect being called the 
"Chiral Magnetic Effect"). This separation of charges 
along the angular momentum vector would be a clear 
T'-violation. 

The expectation from this local T-'-violation is that the 
relative sign of charge separation and angular momentum 
vectors is random in each event. This implies that any 
■p-odd observable should yield zero when averaged over 
many events. An experimental search for this effect must 
therefore involve comparing the measured charge separa- 
tion signal in each event with the expected fluctuations 
due to non-P-violating effects, or equivalently measuring 
correlations among particles in each event. This Letter 
reports the result of such a search performed in 200 GeV 
Au-|-Au and Cu-|-Cu heavy-ion collisions with the STAR 
detector at the Relativistic Heavy Ion Collider (RHIC). 

Observables and Method. In an event, charge sep- 
aration along the angular momentum vector may be de- 
scribed phenomenologically by sine terms in the Fourier 
decomposition of the charged particle azimuthal distri- 
bution 

CX 1 + 2vi COs{cj) - ^lip) + 2V2 cosi2{(j) - irup)) + ... 
+ 2a±sm{(f)~'FRp) + (1) 

where ^I'rp is the azimuthal angle of the reaction plane 
(the plane which contains the impact parameter and 
beam momenta - see Fig.[l]for a schematic), vi and V2 are 
coefficients accounting for the so-called directed and ellip- 
tic flow [5]. The a parameters, a_ — — a_|_, describe the 
"P-violating effect. The spontaneous nature of the pre- 
dicted parity violation means that the sign of and a_ 
vary from event to event and on average (a+) = (i-) ~ 0. 

We may, however, expect non-zero values for the cor- 
relator {aattp) (where a, P represent electric charge -I- or 




FIG. 1: Schematic depiction of the transverse plane in a colli- 
sion of two heavy ions (shown as dotted outlines - one emerg- 
ing from and one going into the page). The azimuthal angles 
of the reaction plane and produced particles with charges a 
and /3 as used in Eqs. [l] and [2] are depicted here. 



— ) since T'-violating contributions to these observables 
would accumulate over many events. One could measure 
(oaO^) by calculating the average (sinZ\(/)Q sinZ\(/)^) over 
all particles of charge a paired with all particles of charge 
(3 (here we have introduced the notation A(j) = (l) — ^pp). 
This is, however, also sensitive to several parity con- 
serving physics backgrounds which produce correlations 
that have nonzero projections along the angular momen- 
tum vector. This led to the proposal 'S^ of the ob- 
servable (cOs((?!)q + 0^ — Ztf'ijp)) = { {cos A(j)a COS Acjip) — 

(sin zi0Q sin Z\0^)) which is sensitive to —{aadp). This 
observable represents the difference between azimuthal 
correlations projected onto the direction of the angu- 
lar momentum vector and correlations projected onto 
the collision reaction plane. By taking this difference, 
these parity conserving correlations only produce back- 
grounds inasmuch as they depend on orientation with 
respect to the reaction plane. As a consequence of the 
(cos Z\0q cos A(t>f)) term, this observable is also sensitive 
to contributions from directed flow and its fluctuations. 
Because the measurements presented in this Letter are 
for a symmetric rapidity region around center-of-mass 
rapidity, these contributions are negligible [TU] . 

In practice, the reaction plane of a heavy-ion collision 
is not known, and one estimates it by measuring the 
"event plane" which is reconstructed from particle az- 
imuthal distributions [5]. In the present analysis, this is 
done by using 3-particle azimuthal correlations in which 
the third particle (labeled "c" below) serves to measure 
the event plane [9l [11] . Assuming that the only corre- 
lation between particles of type "c" and particles a, /3 is 
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through the common correlation to the reaction plane, 
we can evaluate our observable for — {aaap) as 

(cos(0a +(t>P'- 2^Rp)) = (cos(0a + (/)/3 - 20c))/w2,c- (2) 

By varying our choice of the type and momentum range 
of the "c" particles, we have tested this assumption as 
described below. 

Based on available theoretical understanding of the 
Chiral Magnetic Effect we expect the following features 
of the correlator: neglecting any "final state" interac- 
tion with the medium, one expects (a+a+) = (a_a_) = 
— (a+a_) > 0. The likely effect of final state interactions 
in a hot dense medium [5] is a suppression of back-to- 
back correlations (i.e. correlations among two particles 
that travel in opposite directions through the medium), 
causing (a+a_|_) :s> \ (a+a_) |. The dependence of the sig- 
nal on the size of the colliding system has not yet been 
calculated, but one qualitative prediction is that the sup- 
pression of opposite-charge correlations should be smaller 
in collisions of lighter nuclei . For a given collision sys- 
tem, under the assumption that the average size of the 
T'-violating domain does not change with centrality [B], 
the correlator should be inversely proportional to charged 
particle multiplicity, Nch, scaled by a factor accounting 
for the magnetic field in the collisions. Finally, because 
the phenomenon is non-perturbative in nature, we expect 
that the signal should not extend in transverse momen- 
tum far beyond 1 GeV/c, although this may be affected 
by radial flow of the produced particles. 

Data and Detector. The data were obtained with 
the STAR detector [12] during RHIC runs in 2003-04 and 
2004-05. The results are based on 14. 7M Au-|-Au and 
13. 9M Cu-t-Cu events at y/sNN=200 GeV. A minimum 
bias trigger was used with events sorted into centrality 
classes based upon charged particle multiplicity. 

The correlations are reported for charged particle 
tracks measured in the STAR Time Projection Cham- 
ber (TPC) with pseudorapidity \r]\ < 1.0 and transverse 
momentum 0.15 < pt < 2 GeV/c. For event plane de- 
termination, in addition to the main TPC we use mea- 
surements in the two Forward TPCs (2.7 < |?7| < 3.9) 
and two Zero Degree Calorimeter Shower Maximum De- 
tectors (ZDC-SMDs) [131 . The latter are sensitive to the 
directed flow of neutrons in the beam rapidity region. 

The STAR detector is well suited to measure azimuthal 
correlations. The TPC has full azimuthal coverage and 
a charged particle track reconstruction efficiency of ap- 
proximately 85%. Nevertheless, TPC sector boundaries, 
occasional readout channel outages, etc., may introduce 
biases in the analysis. In particular, they may cause in- 
efficiencies that are different for positive and negative 
particles. In evaluating Eq. [2j we correct for detector 
effects (following [HI [15]) by replacing cos(n0i) with 
[cos(n0i) - (cos(n(/)))] for each particle, and similarly for 
the sin(n0i) terms which also appear [TT]. We calculate 
the "re-centering" corrections {cos{n(j))) and (sin(n(/))) as 



a function of time as well as event multiplicity and z- 
vertex position. We also account for the acceptance de- 
pendence on particle ry, pt and charge. Higher order ac- 
ceptance corrections are found to be negligible. 

Experimental Uncertainties. The dominant exper- 
imental systematic error comes from our knowledge of V2 
which is used in Eq. [2] The shaded bands in the figures 
reflect this systematic error with the actual points de- 
termined by applying V2 for TPC particles as measured 
using the reaction plane found in the FTPC. 

Other experimental systematics, including possible bi- 
ases due to acceptance and detector efficiency and errors 
related to track quality cuts, are found to be comparable 
to or smaller than statistical errors [TT]. 

Results. Figure [2] presents {cos{(j)a + 4>p ~ 2^iip)) for 
Au-|-Au and Cu-I-Cu collisions at ^smn=2QQ GeV as 
evaluated using the right-hand side of Eq. [2] (error bars 
indicate statistical errors). The signal in Cu-|-Cu col- 
lisions is larger than the signal in Au-|-Au collisions at 
the same centrality, qualitatively consistent with the ex- 
pected decrease of the signal with increasing multiplic- 
ity. For the Au-|-Au system, opposite-charge correlations 
are clearly smaller in magnitude than same-charge cor- 
relations, in qualitative agreement with the possible sup- 
pression of back-to-back charge correlations. This is sup- 
ported by the observation of a smaller difference in mag- 
nitude between same-charge and opposite-charge correla- 
tions in the smaller Cu-|-Cu system. However, there is a 
large potential background contribution from 3-particle 
clusters to opposite-charge correlations which is discussed 
below and indicated by the thick solid (Au-f Au) and 
dashed (Cu-|-Cu) lines on Fig. [2j 

Figure [3] shows the dependence of the signal on the 
sum of the transverse momentum magnitudes of the two 
particles for the 30-50% centrality range in 200 GeV 
Au-|-Au collisions. We do not observe a saturation or 
drop of the magnitude of the signal at high pt as one 
might naively expect for local V violation. The correla- 
tions are nearly independent of the pt difference over the 
range < \pt,a — Pt,0\ < 2 GeV/c [ll], which excludes 
quantum interference (HBT) or Coulomb effects as pos- 
sible explanations for the signal. We have studied the 
dependence of the signal on \ria — ri/sl [TT], and find that 
the signal has a width of about one unit of rj. 

Physics Backgrounds. We first consider back- 
grounds due to multi-particle correlations (3 or more 
particles) which are not related to the reaction plane. 
This contribution affects the assumption that two parti- 
cle correlations with respect to the reaction plane (L.H.S. 
of Eq. [2]) can be evaluated in practice via three particle 
correlations (R.H.S. of Eq. [2|. Evidence supporting this 
assumption comes from the consistency of same-charge 
results when the reaction plane is found using particles 
"c" detected in the TPC, FTPC, or ZDC-SMD, though 
the FTPC and (particularly) ZDC-SMD analyses have 
large statistical errors in the most peripheral bins. This 
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FIG. 2: {cos{(j)a + 0/3 - 2^Rp)) in Au+Au and Cu+Cu colli- 
sions at \fsNN = 200 GeV calculated using Eq. [5] The thick 
solid (Au+Au) and dashed (Cu+Cu) lines represent HIJING 
calculations of the contributions from 3-particle correlations. 
Shaded bands represent uncertainty from the measurement of 
V2- Collision centrality increases from left to right. 
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FIG. 3: Dependence of (cos(0q, + 0^ — 2tf'ijp)) on ^{pt,a + 
Pt,i3) calculated using no upper cut on particles' pt- Shaded 
bands represent V2 uncertainty. 



multi-particle background should be negligible when the 
ZDC-SMD event plane is used, so it can certainly be re- 
duced and this is an important goal of future high statis- 
tics runs. To study these backgrounds in the current 
analysis, we use the heavy-ion event model HIJING [TS] 
(used with default settings and jet quenching off in all cal- 
culations shown in this Letter) which includes production 
and fragmentation of mini-jets. We find that the contri- 
bution to opposite-charge correlations of three particle 
correlations in HIJING (represented by the thick solid 



and dashed lines in Fig. [2] and |4]) is similar to the mea- 
sured signal in several peripheral bins. We thus cannot 
conclude that there is an opposite-charge signal above 
possible background. The same-charge signal predicted 
by three-particle correlations in HIJING is much smaller 
and of opposite sign compared to that seen in the data. 

Another class of backgrounds (which cannot be re- 
duced by better determination of the reaction plane) con- 
sists of processes in which particles a and P are products 
of a cluster of two or more particles (for example a reso- 
nance decay or jet) and the cluster itself exhibits elliptic 
flow or fragments differently when emitted in-plane com- 
pared to out-of-plane [21 E] . 

For jets with a leading charged particle of pt > 
3 GeV/c, we estimate the contribution using previous 
STAR measurements [T71 [TS] and find it to be negligi- 
ble. To extend the study to lower pt, we rely on HIJING 
calculations of two particle correlations with respect to 
the true reaction plane. These calculations also predict 
the contribution to be small compared to our measured 
signal as shown by the triangles in Fig. |4] 

Resonance decays have the potential to contribute to 
(cos((/)q, + 0/3 — 2^jip)). In addition, previous correlation 
measurements from the ISR [H] and RHIC [201 IH] indi- 
cate that a prominent role in particle production is played 
by clusters. A much smaller signal is expected for same- 
than opposite-charge correlations from resonances, which 
is qualitatively very unlike the signal shown in Fig. [2] 
Kinematic studies demonstrate that it is very difficult 
for the correct sign of fake signal to be created in the 
same-charge correlations without postulating a negative 
value of V2 for the resonances or particles from cluster 
decays. 

To search for other backgrounds to 
(cos((/)q + 0^ — 2tZ/flp)), we have simulated Au+Au 
collisions with heavy-ion event generators MEVSIM [52] , 
UrQMD f237, and HIJING (with and without an elliptic 
flow afterburner implemented as suggested in |8]) for 
comparison and these results (calculated using the 
true reaction plane in all cases) are shown as open 
symbols in Fig. [4] MEVSIM only includes correlations 
due to resonance decays and an overall elliptic flow 
pattern. UrQMD and HIJING are real physics models 
of the collision and so include correlations from many 
different physical processes. Figure [4] shows that no 
generator gives qualitative agreement with data for 
two particle correlations with respect to the reaction 
plane. The models also do not match the measured 
values for reaction plane independent correlations, 

(cos(0a - (j)^)) [n]. 

Other effects we find to produce insignificant contribu- 
tions [TT] include global polarization of hyperons along 
the direction of the system angular momentum. 

Summary. Measurements of three particle correla- 
tions that are directly sensitive to predicted local V- 
violation in heavy-ion collisions have been presented for 
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FIG. 4: (cos(<^c« + 0/3 - 2'1'hp)) results from 200 GeV Au+Au 
collisions are compared to calculations with event gener- 
ators HIJING (with and without an "elliptic flow after- 
burner" ), UrQMD (connected by dashed lines), and MEVSIM. 
Thick lines represent HIJING reaction-plane-independent 
background. 



Au-f Au and Cu+Cu collisions at ^/snn— 200 GeV as a 
function of collision centrality and (in Au+Au) particle 
transverse momentum, with a more complete set of mea- 
surements reported in 

Qualitatively the results agree with the magnitude and 
gross features of the theoretical predictions for local V- 
violation in heavy-ion collisions, but the signal persists to 
higher transverse momentum than expected. The observ- 
able used in our analysis is V-even and might be sensitive 
to non-parity- violating effects. So far, with the system- 
atics checks discussed in this paper and [TT] , we have not 
identified effects that would explain the observed same- 
charge correlations. The observed signal cannot be de- 
scribed by the background models that we have studied 
(HIJING, HIJING+?;2, UrQMD, MEVSIM), which span 
a broad range of hadronic physics. 

A number of future experiments and analyses are nat- 
urally suggested by these results. One of them is the 
dependence of the signal on the collision energy. The 
charge separation effect is expected to depend strongly 
on the formation of a quark-gluon plasma [5], and the 
signal might be greatly suppressed or completely absent 
at an energy below that at which a QGP can be formed. 

Improved theoretical calculations of the expected sig- 
nal and potential physics backgrounds in high energy 
heavy-ion collisions are essential to understand whether 



or not the observed signal is due to local strong parity 
violation, and to further experimental study of this phe- 
nomenon. 
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